Polymer translocation into a fluidic channel through a nanopore 
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Using two dimensional Langevin dynamics simulations, we investigate the dynamics ol polymer 
translocation into a fluidic channel with diameter R through a nanopore under a driving force 
F. Due to the crowding effect induced by the partially translocated monomers, the translocation 
dynamics is significantly altered in comparison to an unconfined environment, namely, we observe 
a nonuniversal dependence of the translocation time r on the chain length A'^. r initially decreases 
rapidly and then saturates with increasing R, and a dependence of the scaling exponent a of r with 
N on the channel width R is observed. The otherwise inverse linear scaling of r with F breaks down 
and we observe a minimum of a as a function of F. These behaviors are interpreted in terms of the 
waiting time of an individual segment passing through the pore during translocation. 

PACS numbers: 87.15.A-, 87.15.H- 



O 



o 
o 



> 

l> 
(N 
(N 

00 

o 

O 



X 



I. INTRODUCTION 

Polymer translocation through a nanopore has at- 
tracted broad interest because it is of fundamental rele- 
vance in polymer physics and is also related to many bi- 
ological processes, such as DNA and RNA translocation 
across nuclear pores, protein transport through mem- 
brane channels, or viruses injecting their DNA into a 
cell. Due to its potentially revolutionary technological 
applications [11,0, including rapid DNA sequencing, gene 
therapy and controlled drug delivery, a considerable num- 
ber of recent experimental [3|-|25| and theoretical [25l - [7l| 
studies have been devoted to this subject. 

The average translocation time r as a function of the 
chain length N is an important measure of the underlying 
dynamics. Standard equilibrium Kramers analysis |72| of 
diffusion across an entropic barrier yields r ~ N'^ for un- 
biased translocation and t ^ N for driven translocation 
(assuming friction to be independent of A^) [l^iUS- How- 
ever, the quadratic scaling behavior for unbiased translo- 
cation cannot be correct for a self-avoiding polymer [42j 
because the translocation time would be shorter than 
the Rouse equilibration time of a self-avoiding polymer, 
TR ^ where the Flory exponent v = 0.588 in 3D 

and i^2D = 0.75 in 2D [zl, |74|- This renders the concept 
of equilibrium entropy and the ensuing entropic barrier 
inap propriate for translocation dynamics. Chuang et al. 
p3 | studied the translocation for both phantom and self- 
avoiding polymers by numerical simulations with Rouse 
dynamics for a 2D lattice model and showed that for 
large N , t ^ N^^'^'^, the same scaling behavior as the 
equilibration time but with a much larger prefactor. This 
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result was recently corroborated by extensive numerical 
simulations based on the Fluctuating Bond (FB) and 
Langevin Dynamics (LD) models with the bead-spring 
approach [49l.l67i. 

For driven translocation, Kantor and Kardar [4^ 
demonstrated that the assumption of equilibrium in poly- 
mer dynamics breaks down even more easily and provided 
a lower bound t ^ N^'^'^ for the translocation time by 
comparison to the unimpeded motion of the polymer. 
Using FB [i^ and LD [43, [5l| models, a crossover from 
T ~ N'^'^ for relatively short polymers to r N^^" for 
longer chains was found in 2D. In 3D, we find that for 
faster translocation processes r ^ N^'^"^ [H, [s^l, while 
it crosses over to r ^ N^'^" for slower translocation, 
such as under weak driving forces and/or of high friction 
constants (55| . Moreover, using linear response theory 
with memory effects and some non-trivial assumptions, 
Vocks et. al. [5^ came up with an alternative estimate 
T N for 3D, which means a = 1.37 in 3D. Their 
a in 3D is consistent with our numerical data for fast 
translocation, but fails to capture the scaling exponent 
for slow translocation. 

However, above physical pictures are based on translo- 
cation into an unconfined trans side. Even for the case 
of an unconfined trans side, the translocated chain is 
highly compressed during the translocation process un- 
der fast translocation conditions jssj , and thus even more 
severe nonequilibrium effects are expected under confine- 
ment. Quite little attention has been paid to the dynam- 
ics of translocation into confined environments. We here 
quantify the effects of the large entropic penalty on the 
confined chain and show that it significantly affects the 
translocation dynamics. In particular we find folded con- 
figurations of the confined chain on the trans side right 
after completion of the translocation process. Moreover 
we show that the force dependence of the translocation 
time significantly changes the 1/F dependence for free 
translocation, and that the almost symmetric form of the 
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FIG. 1: (Color online) A schematic representation of polymer 
translocation through a pore into a 2D confined environment 
under an external driving force F across the pore. The pore 
width is 2a. 

waiting time distribution turns over to a rapid increase 
throughout the translocation process. The paper is or- 
ganized as follows. In section II, we briefly describe our 
model and the simulation technique. In section III, we 
present our results. Finally, we conclude in section IV. 



II. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard-Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction between 
beads is modeled by a short range repulsive LJ potential: 
ULjir) = 4£[(2)i2 _ + £ for r < 2^/^ and for 

r > 2i/V. Here, (T is the diameter of a bead, and e 
is the depth of the potential. The connectivity between 
neighboring beads is modeled FENE spring with 

UpENEir) = -^kRlln{l - r'^/R^), where r is the dis- 
tance between consecutive beads, k is the spring constant 
and Ro is the maximum allowed separation between con- 
nected beads. 

We consider a geometry as shown in Fig. [1] where 
two strips with separation R are formed by stationary 
particles within a distance a from each other. One strip 
( "wall" ) has a pore of diameter 2cr. Between all bead- wall 
particle pairs, there exists the same short range repulsive 
LJ interaction as described above. In the Langevin dy- 
namics simulation, each bead is subjected to conserva- 
tive, frictional, and random forces, respectively, with (76j 
mri = -V{Ulj + Ufene) + Fcxt - ■^v,; + Ff . Here m 
is the bead's mass, ^ is the friction coefficient, v.j is the 
bead's velocity, and Ff is the random force which sat- 
isfies the fluctuation-dissipation theorem. The external 
force is expressed as Fext = Fx^ where F is the external 
force strength exerted on the beads in the pore, and x is 
a unit vector in the direction along the pore axis. 

In the present work, the LJ parameters e, cr, and m 



fix the system energy, length and mass units respectively, 
leading to the corresponding time scale t^j = (mcr^/e)^/^ 
and force scale e/cr, which are of the order of ps and pN, 
respectively. The dimensionless parameters in the model 
are then chosen to be i?o = 2, A: = 7, ^ = 0.7, and 
F = 0.5...15. 

In our model, each bead corresponds to a Kuhn length 
(twice of the persistence length) of a polymer. For a 
single-stranded DNA (ssDNA), the persistence length of 
the ssDNA is sequence and solvent dependent and varies 
in a wide range, to our knowledge, usually from about 
1 to 4 nm. We assume the value of a ~ 2.8 nm for a 
ssDNA containing approximately four nucleotide bases. 
The average mass of a base in DNA is about 312 amu, so 
the bead mass m k, 1248 amu. We set fc^T = 1.2e, which 
means that the interaction strength e is 3.39 x 10~^^ J at 
actual temperature 295 K. This leads to a time scale of 
69.2 ps and a force scale of 1.2 pN. Each base (nucleotide) 
is estimated to have an effective charge of 0.094e from 
Ref. 0, leading to an effective charge of a bead with 
four bases of 0.376e. Thus, the voltage across the pore is 
between 28.1 and 843 mV for varying F from 0.5 to 15, 
within the range of experimental parameters [H-IH . 

The Langevin equation is then integrated in time by a 
method described by Ermak and Buckholz [77| . Initially, 
the first monomer of the chain is placed in the entrance 
of the pore, while the remaining monomers are under 
thermal collisions described by the Langevin thermostat 
to obtain an equilibrium configuration. Typically, we 
average our data over 1000 independent runs. 

In Fig. [1] we assume that both ends of the channel are 
open. If that were not the case, such as for virus capsids, 
on the entering of the chain the fluid in the conflned 
region would necessarily have to leave. To address such 
a problem it is therefore relevant to explicitly include the 
solvent in the model, which is beyond the scope of the 
present study. 



III. RESULTS AND DISCUSSION 

Consider a polymer, such as DNA, confined to a 
nanochannel of width R with R being less than the radius 
of gyration of the molecule. The response of the poly- 
mer to confinement is primarily dictated by the relative 
value of R with respect to the chain persistence length. 
Depending on whether R is larger (de Gennes regime) or 
smaller (Odijk regime) than the chain persistence length, 
different scaling behaviors of the longitudinal size of the 
chain, R\\ , as a function of R were predicted in the pio- 
neering theoretical studies by dc Gennes [t^ and Odijk 
[TSj . respectively. In the Odijk regime, the physics is 
dominated not by excluded volume but by the interplay 
of confinement and intrinsic polymer elasticity. In this 
work, we only consider the de Gennes regime, where the 
blob picture [73| is valid. To consider the Odijk regime, 
we would need to take into account the chain stiffness 
in the model. While this is certainly interesting we here 



3 




FIG. 2: Translocation time r as a function of F for 7? = 3.5 
and R = oom 2D. The chain length is = 128. 



focus on the polymeric aspects in the flexible chain limit. 

According to the blob picture, for a polymer confined 
between two strips embedded in 2D the chain will ex- 
tend along the channel forming blobs of size R. Each 
blob contains g ~ [Rj a)^!^'^" monomers, where is 
the Flory exponent in 2D [tI, [tI], and the number of 
blobs is Uh = N/g = N{a / Rf/"^'' . Then, the blob 
picture predicts the longitudinal size of the chain to be 
~ rifei? ~ Nai^f/"^"-^ - iVi?-i/3 [7l,[7|. The lon- 
gitudinal relaxation time T|| is defined as the time needed 
by a polymer to move a distance of the order of its longi- 

tudinal size, R\\ . Thus, T|| scales as T|| " 

with D = 1/N being the diffusion constant. The free 
energy cost in units of ksT is J" = Nia/RY^"^" . 

For polymer translocation into confined environments, 
a driving force is necessary to overcome the entropic re- 
pulsion f{R) exerted by already translocated monomers. 
Due to the highly non-equilibrium property of the 
translocation process, it is difficult to estimate the resist- 
ing force f{R). If we assume that, for slow translocation 
processes, the resisting force f{R) scales as /(i?) = CR'^ , 
with C and 7 being the associated prefactor and the 
scaling exponent, respectively; then, under an external 
driving force F in the pore the translocation time r can 
be written as r ^ p^fi^^j ^ f^i-cr-i /f) '^ith a being 
the scaling exponent of r with chain length A^. Due to 
Too ~ ^ for an unconfined system with R 
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have 1 — ^ ^ CR? jF. Based on this relationship, we 
examine the dependence of r on R. 



Translocation time as a function of the driving 
force 



As shown in Fig. [5] for unconfined system R = 00, 
the dependence of r on the driving force scales as F~^ 




FIG. 3: (Color online) Waiting time distribution for A'' = 
128 and R = 3.5 under different F. The waiting time of 
monomer s is defined as the average time between the events 
that monomer s and monomer s + 1 exit the pore. 




FIG. 4: (Color online) Typical chain conformation at the mo- 
ment just after translocation for A'' — 128 and R = 3.5 under 
different F. 



for weak driving forces F < 4. This simple scaling be- 
havior can be understood by considering the steady state 
motion of the polymer through the nanopore. The av- 
erage velocity is determined by balancing the frictional 
damping force (proportional to the velocity) with the ex- 
ternal driving force. This leads to an average velocity 
proportional to the driving force F, and hence a translo- 
cation time T ^ F^^. For F > 4, r ^ 
pronounced non-equilibrium situation where the chain is 
highly distorted, as illustrated in our previous work jss} . 
For this case, only part of the chain on the cis side can 
respond immediately, while the remaining part near the 
chain end does not feel the force yet. As a result, a part 
of the chain on the cis side is deformed to a trumpet and 
even stem-and-flower shape, while the translocated por- 
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tion on the trans side has a compact spherical shape, as 
it does not have time to diffuse away from the pore exit. 
In a recent theoretical study [H^ , the effect of a trumpet 
shape of the chain on the cis side, on the translocation 
dynamics was found to cause a breakdown of the t ^ F^^ 
scaling. However, this theory neglects effects due to the 
compacted chain structure on the trans side. 

However, for confined systems the effect of the driv- 
ing force on the translocation time is completely differ- 
ent, as shown in Fig. [2] for i? = 3.5. With increas- 
ing F, the translocation time initially decreases rapidly, 
and then almost saturates for F > 10. As noted above, 
F{i-CRri /F) equilibrium process. With de- 

creasing F, the resisting force /(i?) is more and more 
important and it slows down translocation. This is the 
reason why the initial slope is faster than that for translo- 
cation into a free trans side. With increasing F, f{R) 
becomes important again, which indicates that the re- 
sisting force f{R) induced by crowding effects from al- 
ready translocated monomers plays an important role in 
the observed behavior. 

To further understand this behavior, we examine the 
dynamics of a single segment passing through the pore 
during translocation. The nonequilibrium nature of 
translocation has a significant effect on it. We have nu- 
merically calculated the waiting times for all monomers 
in a chain of length N. We define the waiting time of 
monomer s as the average time between the events that 
monomer s and monomer s -I- 1 exit the pore. In our 
previous work, we found that the waiting time depends 
strongly on the monomer positions in the chain under 
the driving force F = 5 [H, [4^. For relatively short 
polymers, such as N = 100, the monomers in the mid- 
dle of the polymer need the longest time to translocate 
and the distribution is close to symmetric. In contrast, 
for the confined system with i? = 3.5 and N ~ 128, we 
find a slow increase after the monomers in the middle 
of the polymer (i.e., beyond s = N/2) and a further, 
more rapid increase after s « 100 for F = 5, see Fig. 
[3l Increasing the driving force further to F = 15, the 
waiting time rapidly increases throughout. Moreover, it 
takes longer for monomers at s > 100 to exit the pore 
for strong driving force = 15 than that for F = 5. 
These results indicate that during later stages of translo- 
cation the high density of segments in the channel slows 
down the translocation and the inverse proportionality 
T ~ F~^ breaks down. Especially for strong driving 
forces the crowding effect is more pronounced, leading 
to almost the same translocation time for = 10 and 
F — 15. These results are different from translocation 
into two parallel walls (3D), where we observed a rela- 
tively fast turnover between the scaling r ^ ^-i-25 
weak forces and the behavior r ~ F^^.s strong forces 

Fig. |4] shows the typical chain conformation at the mo- 
ment just after translocation for = 128 and R = 3.5 
under different F. Compared with the equilibrated chain, 
the chains become more compressed with increasing F. 
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FIG. 5: Translocation time t as a function of the chain length 
A'^ for R = 3.5 and different F in 2D. The insert shows the 
scaling exponent a of r with A for different driving forces 
F. For F=0.5, 3, 5, 10 and 15, a = 2.54 ± 0.04, 1.85 ± 0.01, 
1.88 ± 0.01, 2.03 ± 0.02, and 2.21 ± 0.04, respectively. 



In particular, a folding of the chain is observed for F > 2, 
reflecting stronger resisting forces with increasing F . At 
_F > 10 there even occur triple layers in the chain configu- 
ration. As can be seen the chains in this folded state are 
almost entirely devoid of the larger-amplitude undula- 
tions observed in the equilibrium configuration. This also 
provides a way to induce chain folding by driving poly- 
mers into unidimensionally confined environments such 
as small channels. 

We stop to note that here we confine the analysis to 
the variation of the driving force F . As demonstrated in 
Ref. [m the variation of the friction coefficient ^ has a 
similar effect. 



B. Translocation time as a function of the chain 
length 

Previously [H, |4^ [5l|, we investigated the polymer 
translocation into an unconfined trans side (R = oo) 
under an external driving force = 5 in the pore. 
Both the 2D fluctuating bond model with single Monte 
Carlo moves [l^ and Langevin [i^ [5l| simulations show 
r ~ N^" ^ N^-^ for relatively short chains N < 200. 

However, for polymer translocation into the region be- 
tween two strips, the dynamics is completely different. 
As shown in Fig. [Sj for the same driving force F — 5 
for R = 3.5, we flnd r - jvi-^sio.oi transloca- 
tion velocity v ^ jYi.08±o.02 ^^j^Qgg i^j^g inversely lin- 
ear scaling v ^ N^^. During the translocation process, 
the chain moves a distance of Rg + instead of 2Rg 
for an unconfined system, where is the longitudi- 
nal size of the chain. For a polymer confined between 
two strips embedded in 2D, the longitudinal size of the 
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FIG. 6: (Color online) The radius of gyration of the chain 
before translocation and at the moment just after the translo- 
cation for R = 3.5 and different F. Here, the x direction is 
perpendicular to the wall, and y is the direction along the 
wall. 



chain is _R|| ^ N (see above scaling results from the blob 
picture), which is the dominant term over the second 
term, Rg ^ N'^^° . Thus, the translocation time roughly 

scales as T ^ N'^. Decreasing the driving force to 

F = 3.0, the scaling exponent almost does not change, 
T ~ jyi.85±o.oi Decreasing the driving force further 
to = 0.5, we find the translocation dynamics crosses 
over to another regime with a = 2.54 ± 0.04 at least for 
N < 64, which is still lower than the expected scaling 
exponent 3.0 of the relaxation time as a function of the 
chain length for a chain in a narrow channel. Due to ex- 
pensive computation, we cannot access cases with longer 
chains and/or weaker driving forces. Increasing the driv- 
ing force F from 5 to 10 and 15, we find r ^ jy^ °^='=°"^ 
and T ^ j\r2.2i±o.o4^ respectively. This is, however, an a 
priori unexpected increase of a with increasing F, which 
indicates that the translocation is slowed down for strong 
driving forces. Particularly, for N = 256 the transloca- 
tion time almost does not change for F ^ 5, 10 and 15, 
due to the formation of a folded chain conformation. 

Altogether, compared with an unconfined environment 
we observe the nonuniversal dependence of the transloca- 
tion time T on the chain length N. Particularly, there is 
a minimum of a as a function of F, as shown in the insert 
of the Fig. [51 due to the more severe crowding effect for 
stronger driving forces. The scaling exponent does not 
change for _F = 15 compared with F = 5. These results 
are different from the translocation into two parallel walls 
(3D) with the same R, where for weak driving forces r 
scales exactly in the same manner as the chain relaxation 
time and crosses over to exponent 1.37 for strong driving 
forces [ill. As intutitively expected, these observations 
demonstrate that confinement effects are more relevant in 
the ID channel configuration compared to the 2D con- 
finement of the chain sandwiched between two parallel 
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FIG. 7: Waiting time distribution for different A'' under the 
driving force F = 5 and R = 3.5. 



walls. 

Fig. [5] shows the radius of gyration of the chain be- 
fore the translocation and at the moment just after the 
translocation for R = 3.5 and different F. Just after 
the translocation, the chains are compressed in the y 
direction compared with the equilibrated chain. How- 
ever, the scaling exponent of Rg^y ^ N almost docs not 
change even for F = 15. Fig. [7] shows the waiting time 
distribution for different chain lengths for F — 5.0 and 
R — 3.5. Similar distributions are observed for different 
chain length. This behavior reflects the crowding effect 
of partially translocated monomers, which greatly slows 
down translocation. 



C. Translocation time as a function of the channel 
width R 

As shown in Fig. [51 there exist two regimes in the be- 
havior of T as a function of R for F = 5.0 and N = 128. 
T decreases rapidly with increasing R and then almost 
saturates for larger R. The insert of Fig. [51 shows 
1 — — as a function of R. For i? < 6, we find the 
exponent 7 = —1.06 ± 0.07. For equilibrium transloca- 
tion process f{R) which gives the exponent 
~^/v2D = —1.33. The numerical exponent 1.06 does not 
agree with this result, which indicates that the translo- 
cation is a highly non-equilibrium process. 

Fig. [9] shows the radius of gyration of the chain before 
translocation and at the moment just after the translo- 
cation for N = 128 and F = 5 versus the channel size 
R. Just after the translocation, the chains are highly 
compressed along the channel and Rg^y ^ R"^-^'^ for 
R < 5 compared with the equilibrated chains, where 
Rg^y ^ R^^-^"^ as predicted. 

Fig. [TUl depicts the waiting time distribution for differ- 
ent R. For larger R, wc observe a symmetric distribution 
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FIG. 8: Translocation time r as a function of R for chain 
length A'' — 128 under the driving force F — 5. 
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FIG. 10: Waiting time distribution for different R. The chain 
length A*' = 128 and the driving force F = 5. 
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FIG. 9: (Color online) The radius of gyration of the chain be- 
fore translocation and at the moment just after the transloca- 
tion for TV = 128, F = 5 versus the channel size R. Here, the 
X direction is perpendicular to the wall, and y is the direction 
along the wall. 



with respect to the middle monomer s = N/ 2 [4^ |4^ . 
With decreasing R, the waiting times increase. Particu- 
larly, it takes much longer time for monomers s > N/2 
to exit the pore. For i? = 5, the waiting times approxi- 
mately saturate after s > N/2, while they continuously 
increase for i? < 4. This behavior is due to the crowding 
of the translocated monomers. As to the first moment 
of the translocation coordinate (s(t)) ^ , previous re- 
sults [H show /3 sa 1/a = 0.67 with a « = 1.50 for 
R = oo. Fig. [n] shows {s{t)) - fl^ for N = 128 and 
F = 5 for different R. For i? = oo, ^ = 0.70 ± 0.01 as 
expected and it continuously decreases to 0.58 ± 0.01 for 
i? = 3. Here, /3 < 1 is a signature of anomalous diffusion 
[t^. The values of /3 in Fig. [TT] imply that the scaling 
exponent a can change from « 2i/ to « 2, depending on 



100 



V 



10 



'1 1 1 ' ' 

F=5 






N=128 






Slope: 0.70+0.01 








— □ — 


R=oo 




— o — 


-R=9 : 


Slope: 0.58±0.01 


— A — 


R=7 : 




— O — 


R=6 : 




— o- 


R=5 . 




— <] - 


R=4 




— [>- 


R=3 - 


'1 1 1 ' ' 

1 10 100 

time 




1000 



FIG. 11: The number of translocated monomers as a function 
of time, {s{t)), for different R. The chain length A'^ = 128 and 
the driving force F = 5. 
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FIG. 12: Translocation time as a function of the chain length 
for different R. The driving force F = 5. 
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R due to highly non-equihbrium effects. This is clearly 
indicated in Fig. [H where a = 1.52 ± 0.02, 1.79 ± 0.02 
and 1.93 ± 0.02 for i? = cxo, 5 and 3, respectively. 

IV. CONCLUSIONS 

Using two-dimensional Langevin dynamics simula- 
tions, we investigated the dynamics of polymer translo- 
cation into a narrow channel of width R through a 
nanoporc under a driving force F. Due to the crowding 
effect induced by the partially translocated monomers, 
the translocation dynamics is greatly changed compared 
with an unconfined environment. Namely we observe a 
nonuniversal dependence of the translocation time r on 
the chain length N: r initially decreases rapidly and then 
saturates with increasing R, and, moreover, an R depen- 
dence of the scaling exponent a in the law t ^ is 
observed. The inversely linear scaling of r with F breaks 
down and we observe a minimum of a as a function of 
F. These behaviors are clearly interpreted in terms of 
the waiting time of an individual segment for passing 
through the pore during translocation, as well as in terms 
of statistical quantities such as the width of the chain. In 
particular we observe folded configurations of the chain 



in the channel, at high driving forces even a structure 
with a double-fold exists. Similar effects of nonequilib- 
rium nature are expected in biological cells. There, the 
passage of a translocating chain is opposed by a crowded 
environment in which larger biomolecules occupy more 
than 40% of the volume [sO, IHi. 

Our findings are of interest from a purely polymer 
physics point of view, contributing to a more com- 
plete picture of polymer translocation through a narrow 
pore, once more pronouncing the importance to consider 
nonequilibrium situations. A direct applicatoin of our 
findings may be of relevance for parallel fluidic channel 
setups consisting of stacked parallel slices. In each slice 
the chain initially is in an unconfined two-dimensional 
conformation, before being forced through the pore into 
the effectively ID channel. 
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